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Contrary to a common argument that a small increase in the strength of the strong force would 
lead to destruction of all hydrogen in the big bang due to binding of the diproton and the dineutron 
with a catastrophic impact on life as we know it, we show that provided the increase in strong 
force coupling constant is less than about 50% substantial amounts of hydrogen remain. The reason 
is that an increase in strong force strength leads to tighter binding of the deuteron, permitting 
nucleosynthesis to occur earlier in the big bang at higher temperature than in the standard big 
bang. Photodestruction of the less tightly bound diproton and dineutron delays their production 
to after the bulk of nucleosynthesis is complete. The decay of the diproton can, however, lead to 
relatively large abundances of deuterium. 

PACS numbers: 26.35.+c,98.80.Bp 
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I. INTRODUCTION 

The weak anthropic principle has often been used to 
infer hmits on the range of the strong force strength or 
coupling parameters consistent with life as we know it. A 
common argument is that a small increase of the strong 
force strength will bind the dineutron and the diproton, 
leading to a large increase in the rate of the p+p and 
n+n reactions, so that big-bang nucleosynthesis (BBN) 
leads to all the protons being converted into isotopes of 
helium, leaving no hydrogen necessary for chemistry vital 
to life. For example, Dyson ^ writes "If a helium- 2 nu- 
cleus could exist, the proton proton reaction would yield 
a helium-2 nucleus plus a photon, and the helium-2 nu- 
cleus would in turn spontaneously decay into a deuteron, 
a positron and a neutrino. As a consequence there would 
be no weak interaction hang-up, and essentially all of the 
hydrogen existing in the universe would have been burned 
to helium even before the first galaxies had started to 
condense." Barrow and Tipler [3| state on p322 of their 
book "If the strong interaction were a little stronger the 
diproton would be a stable bound state with catastrophic 
consequences all the hydrogen in the Universe would have 
been burnt to IIe2 during the early stages of the Big 
Bang and no hydrogen compounds or long-lived stable 
stars would exist today. If the diproton existed we would 
not!" 

In this paper, we consider how BBN is altered by the 
existence of bound diproton and dineutron nuclei, tak- 
ing into account some physical processes that so far have 
been overlooked. To relate the binding energies of the 
diproton, the dineutron and the deuteron to the relative 
strength of strong force, we use the same square well 
potential model as Barrow [3], who found that a 9% in- 
crease in the strong force coupling constant, as, is suffi- 
cient to bind the dineutron and a 13% increase will bind 



the diproton. The needed increase in coupling constant 
to bind the diproton was confirmed by Pochet et al [1] for 
a more realistic nuclear potential. For small increases in 
as, the diproton and dineutron binding energies are suffi- 
ciently small that photodestruction will prevent buildup 
of large amounts of these isotopes before freeze out oc- 
curs. 

In section 2, we briefly review the physics of BBN rel- 
evant to our investigation of the effects of increased as ■ 
In section 3 we describe how we determine the needed 
rates for reactions involving dineutrons and diprotons. 
Due to the electrostatic repulsion of the protons, there 
is a narrow range of strong force strength for which the 
dineutron is bound but the diproton is not. Nucleosyn- 
thesis in this regime is considered in section 4. In section 
5, we consider the regime in which the diproton is also 
bound. Finally, in section 6 we give our conclusions. 



II. BBN IN THE RADIATION DOMINATED 
ERA 

In the standard hot big bang Q , nucleosynthesis takes 
place in the radiation dominated era that begins after 
electron pair annihilation is complete. At the begin- 
ning of this era the radiation and matter temperature 
is ^ 4 10^ K. Although the p + n d + 7 reaction is 
rapid, due to the relatively small binding energy of the 
deuteron, photodestruction prevents significant amounts 
of being formed until the temperature has dropped 
to ~ 10^ K. Further reactions continue, building up an 
appreciable amount of ''He and traces of other light ele- 
ments, until freeze out occurs when the Universe is about 
15 minutes old and has a temperature of ~ 5 10* K. Dur- 
ing the radiation era, the temperature varies with time, 
t. as 



Tg = 13.8 f 



-1/2 



(1) 



• ^j immacd@udel . edu| 



where Tg is the temperature in units of 10^ K and t is 
measured in s. The baryonic density in g cm~3 is ap- 
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proximately given by 



p = 3.3 10\Tg^ 



(2) 



where 77 is the ratio of the number of baryons to the 
number of photons. To explore the effects of bound 
diproton and dineutron on BBN, we have written a 
computer program that includes all the important 
reactions [6] 7] plus a few more relevant to universes in 
which the diproton and dineutron are stable. Because 
the binding energy of the deuteron is increased when 
Us is increased, nucleosynthesis can begin at higher 
temperatures than for standard BBN. We have modified 
equations (1) and (2) to include conditions in the 
leptonic era. We take Tg — 50 as the initial temperature. 
The initial abundances of neutrons and protons are set 
to their equilibrium values, with all other abundances 
set to zero. The value of -q is taken to be 4 10"^", 
which gives for the standard big bang good agreement 
with observed light element abundances. As a test of 
our code, we have compared our standard big bang 
nucleosynthesis results with those obtained with the 
public_bigbang code (which can be downloaded from 
,http: / / cococubed.asu.edu/code_pages/net-bigbang . shtml ) . 



III. ADOPTED RATES FOR REACTIONS FOR 
REACTIONS INVOLVING DINEUTRONS AND 
DIPROTONS 

The production of diprotons and dineutrons in the big 
bang will depend on the competition for neutrons and 
protons between the p-\-p ^ pp and n + rt — > nn + 7 
reactions and the p + n ^ d + 7 reaction. To determine 
the rates of the p + p ^ + 7 and n + n — > nn + 7 
reactions, we make use of known results for the p + p ^ 
d + + v reaction. We assume that the intrinsically 
nuclear part of the interaction potential is the same for 
all three reactions. To account for the difference in rates 



of the p+p — > d+e^+u reaction in which the intermediate 
state experiences a weak decay and the two reactions in 
which the intermediate state decays electromagnetically, 
we introduce a factor fw-e- Using the p + p S-factor @], 
a straightforward integration for the n + n — > nn + 7 
reaction gives for the reaction rate per particle in cm^ 
mol~^ s~^ 

Na{ctv) = /^_el0-''(1.78 + 1.80r9+1.99r2)/Tg'/2 (3) 



For the rate of the p + p — 
multiply the p + p — > d + 



pp + J reaction, we simply 
+ u rate by /^,_e, 



NA{crv) = /„_el0"^^(4.08 + I5.6T9 + e.lGTg^ + 0.588T| 
-0.0465Tg^)e-3-3«i/^«''VT9/' (4) 

At first sight it might be thought that fw-e could sim- 
ply be obtained from consideration of the p + n — s- d + 7 
reaction. However the comparison is complicated by the 
effects of nucleon spin on the nuclear interaction. The 
stability of the deuteron is a result of the nuclear force 
being stronger when the nucleons have parallel spin than 
when the spins are opposite. The Pauli exclusion princi- 
ple requires that the nucleons in the diproton and dineu- 
tron have opposite spin. We can roughly estimate fw-e 
by considering the nuclear S-factors for similar reactions. 
For example, comparing 5*0 = 2.5 10~^ KeV barn for 
d(p,7)3He with 5*0 = 3.8 IQ-^^ KeV barn for p(p, P+v)d 
,9j indicates that fw-^e ^ 10^^. The nuclei in this re- 
action do not have the same spins as the nuclei in the 
reactions of interest. However a similar result is ob- 
tained by considering a reaction with the correct spins, 
^^N-|-p — >^'*0-)-7. The non-resonant contribution to this 
reaction has 5*0 ~ 3 10~^ KeV barn[l|. In figure 1, we 
plot the ratios (CTw)„p/(crw)„„ and {(tv)^^/ {av)^^ against 
Tg for the case in which a value of 10^* is adopted for 
fw-e- For the p + n ^ d + ^ reaction rate, we have used 



J 



Na{(tv) = 4.42 10^(1.0 + 3.75Tg + imT^ + 0.747T| + 0.0197r9^ -t- 3.00491 10"^r|)/ 



(1.0 + 5.47T9 + 5.62rg2 



0.489r93 



7.47 lO^^Tg') 



r 



(5) 



We see that at temperatures relevant to BBN, dineutron 
and diproton production will not be important unless 
fw-e ^ 10^^. Because of the uncertainties involved in 
estimating fw-ei in the following we consider a range of 
fw-e values. 

The dineutron and diproton production will also de- 
pend on the rates of the reverse reactions pp + 7 p + p 
and nn+7 n+n. For small increases in the strength of 
the strong force coupling constant, the binding energies 
of dineutron and diproton will be of order kT during the 



nucleosynthesis phase of the big bang. Hence, the thresh- 
old value of e^T'/'^*, where is the photon energy, is not 
large compared to unity, and the usual approximation 
that the reverse rate is the forward rate multiplied by a 
factor 



A. 



= 7.07 10^ t!^^ e-Q/'=^ 



(6) 



NA{av) 

cannot be used indiscriminately. Although the rate of the 
nn + "f n + n reaction can be evaluated analytically in 
terms of Debye functions, we find it simpler to evaluate 
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FIG. 1: Dependence of the ratios {'^'v)„p/{o'v}^^ and 
i'^'") npK'^'") pp temperature in units of 10^ K for fw-e, = 

this rate numerically, together with the + 7 p + p 
reaction rate. For a binding energy Q — 150 keV, we find 
that the above approximation underestimates the reverse 
rate by only 13% at Tg = 4. We determine a correction 
to the approximate rate by multiplying it by a factor of 
form 1/(1 + aiq + 02(7^ + azq^) where q = e~'3/fcT^ r^j^g 
coefficients ai, 02, 03 are found by fitting to the numerical 
results. An important consequence of an increase in 
is that the binding energy of the deuteron will also be 
increased. This will reduce its photodestruction rate and 
allow production to occur at higher temperature than 
in the standard big bang. We take this increase in the 
deuteron binding energy, Qd into account in calculating 
the photodestruction rate of ^H. 

We determine the binding energies, Qd, Qnn and Qpp 
by using a square well potential for the nucleon inter- 
action of depth V and radius r. To characterize the 
strength of the strong force, we use the relative strong 
charge, G, which is related to the strong force coupling 
constant by 

where go is the standard strong charge value. The depth 
of the potential is then 

V = G^Vo (8) 

We use the same values for Vq and r as Barrow Q. For 
the deuteron J'^ = 1+ ground state Vb = 36.2 MeV and 
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FIG. 2: Dependence of the deuteron, dineutron and diproton 
binding energies in MeV on the relative strong charge G. 

r — 2.02 fm. For the J'^ = 0+ ■^H and dineutron states, 
Vo = 14.0 MeV and r = 2.59 fm. The diproton poten- 
tial is taken to be 0.56 MeV shallower than the dineutron 
potential to account for the Coulomb energy. The depen- 
dences of the resulting binding energies on G are shown 
in figure 2. 

We now consider the rates of Icptonic transformations 
between dineutrons, diprotons and deuterons. For small 
binding energies, we expect the time for dineutron de- 
cay to be comparable to that of the neutron. Similarly, 
since the proton is stable, the lifetime of the diproton is 
likely to be quite long. A comparison with the leptonic 
rates between n and p [ll| indicate that the enhance- 
ments of the overall rates, including electron capture, at 
BBN temperatures will be modest and hence for simplic- 
ity we neglect these enhancements. To estimate the weak 
rates, we use the ft factors for the corresponding tran- 
sitions in the beta decays of the analog nuclei ^^O and 

The major ^''Q decay channel is to an excited J'^ = 0+ 
state 2.313 Mev above the ground state of i**N (99.3%, 
log ft = 3.4825). This is analogous to diproton decay to 
the spin singlet state of ^H. There are also decays to the 
l"*" ground state (log ft = 7.279) and to a 1"*" excited state 
at 3.948 MeV (log ft = 3.131). We have calculated the 
Fermi integrals with the relativistic form of the Fermi 
factor for the two higher energy decay channels which 
correspond to those of the diproton. For the diproton 
decay channel corresponding to the dominant ^^O decay, 
we find that an accurate approximation to the decay time 
scale is 

log ti/2 = 3.490 - 5 log E+ (9) 
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where the maximum positron kinetic energy measured in 
MeV is, in terms of binding energies, 

E+ ^ Q,m-Qpp- (TO,i-TOp)-771e = Q„„ - Qpp - 1 .8043. 

(10) 

For the decay to the ground state 

logti/2 = 7.287-51ogS+ (11) 

where now 

E+ ^ Qd~ Qpp - [mn - m.p) - nrie ^ Qd- Qpp - 1.8043. 

(12) 

The analog of the dineutron, ^^C, decays only to the 
ground state of ^^^N with half-hfe ii/2 = 1.8 lO" s (log ft 
— 9.040). We find that the half-life for the correspond- 
ing dineutron decay is approximately related to the max- 
imum electron kinetic energy in MeV by 

= (63/i;_)5 (13) 

where 

= Qd- Qnn - {rUn -nip-me) ^ Qd^ Qnn ^ 0.7823. 

(14) 

Since we expect the dineutron will also have a decay chan- 
nel to the 0+ excited state of (which will exist if the 
dineutron is bound), we use for this channel the ft fac- 
tor for the corresponding transition in ^^O. The energy 
difference between the dineutron and the singlet state 
will always be about the difference in mass of the neu- 
tron and proton, 1.3 MeV, and the /-factor will then be 
about 1.8, which gives ti/2 = 1.5 10'^ s. Combining the 
rates for the two decay channels gives 

A„„ ^ 4.62 10"'' + {E_/58.5f. (15) 

The resulting decay time scales, A~^, are plotted against 
G in figure 3. 

Finally, we need to also consider additional reactions 
that arise when the diproton and dineutron are bound. 
The most rapid reactions are likely to he pp + n d + p 
and nn+p —>■ d + n. Due to the complexity of calculating 
reaction rates even for few nucleon systems (see for ex- 
ample Marcucci et al. [13 )i we settle for estimating when 
these two reactions are likely to be important by compar- 
ing the results of calculations in which these reactions are 
completely neglected with the results of calculations in 
which they are assumed to be instantaneous. 

IV. THE n+n nn REACTION REGIME 

Due to the electrostatic repulsion of the protons, there 
is a narrow range of strong force strength for which the 
dineutron is bound but the diproton is not. According 
to the square well potential model, this range is 1.043 < 
G < 1.063. The dineutron binding energy in this regime 
is Qnn = 0—14 KeV. The relevant new reactions for BBN 
are n + n — > nn + j, nn + j — > n + n and nn d + e~ +v. 
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FIG. 3: Dependence of the dineutron and diproton beta decay 
life times in seconds on the relative strong charge G. 

In figure 4, we show for G = 1.06 how the final mass 
fractions of ^H, ^He and the dineutron depend on fw~e 
when only the n + n nn + 7 reaction is included. We 
see that there is no significant production of dineutrons 
unless fw-e ^ 10^^. Also the H abundance increases with 
fw-e because the n + n —f nn + 7 reaction removes the 
neutrons before the n+p ^ d + e'^ + 1^ reaction can take 
place. 

In figure 5, we show the final dineutron abundance 
when the nn -f 7 — > n -f n reaction isalso included for dif- 
ferent values of G. Dineutron production is small unless 
its binding energy is comparable to that of the deuteron. 
Otherwise, the neutrons are removed by proton capture 
before significant amounts of dineutron can be produced. 
For values of G large enough to give the dineutron a 
binding energy greater than 2 MeV, the diproton would 
certainly be bound. 

The n -|- n — > nn + 7 reaction could lead to significant 
production of if the dineutron leptonic decay occurs 
more quickly than its photodestruction. In figure 6, we 
show how the final mass fraction of 'H depends on the 
time scale of the nn d+e~+v reaction for G = 1.06. In 
order to have significant dineutron production to occur, 
we have set fw-e = 10^°. We see that, if fast enough, the 
leptonic decay increases the hydrogen abundance. This 
is because the set of reactions 

n + n ^ nn + 7 
nn —fd + e~+v 
d + ^ n + p 
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FIG. 4: Final mass fractions of (solid line), *He (long dash 
line) and dineutron (short dash line) when only the n+n 
nn reaction is included. 








-28 ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' 
5 7 9 11 13 15 17 19 21 23 25 

log t a 

FIG. 5: Final abundance of dineutron for G values from bot- 
tom to top of 1.25, 1.50, 1.75 and 2.00. 

converts a neutron into a proton. For this to happen, 
the dineutron decay must occur on a time scale of 10^^^ 
s or less, which is much less than the estimate above, 
Tnn ~ 10^ s. Hence it is unlikely that this set of reactions 
is important. 

Finally we consider the effects of the nn + p ^ d + n 
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FIG. 6: Dependence of hydrogen final abundance for G = 1.06 
on dineutron life time (in s). 

reaction. To gauge the importance of this reaction, we 
assume that it is instantaneous. The final mass fractions 
of ^H, "^He are shown in figure 7, again for G = 1.06. We 
see that this reaction leads to small reductions in the H 
abundance for fw-e > 10^'^. Hence this reaction does not 
have a major effect on BBN when only the dineutron is 
bound. 

To summarize the results presented in this section, we 
find that for values of the strong force coupling constant 
at which the dineutron is bound and the diproton is un- 
bound, there are no catastrophic impacts on BBN. 



V. THE p-hp ^ pp REGIME 

For G > 1.065 both the diproton and the dineutron 

are bound. Increased G also binds the deuteron more 
tightly, allowing it to be formed earlier in the big bang at 
higher temperatures, where the less tightly bound dipro- 
ton and dineutron are easily destroyed by energetic pho- 
tons. We first consider only the effects of increased G on 
deuteron binding by setting fw-e = 0. Figure 8 shows 
how the final hydrogen and helium abundances depend 
on G. In the standard big bang ^H production begins 
when the temperature has dropped to about 10^ K. For 
G > 1.2, Qa is high enough that ^H production begins 
in the leptonic era. The final ^H abundance is then ap- 
proximately the difference in the equilibrium proton and 
neutron abundances at the temperature at which pho- 
todestruction of ^H becomes unimportant. 

When fw-e > 0, if diproton production occurs it does 
so long after the primordial neutrons have been consumed 
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FIG. 7: and ''He mass fractions when the nn + p ^ d + n 
reaction is included. 
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FIG. 9: Final mass fractions of n, ^H, ^H, ^He and '*He for 



G = 1.3 and 



10 to illustrate typical behavior. 
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FIG. 8: Dependence of the final abundances of H (solid line) 
and ''He (broken line) on G in the absence of production of 
dineutrons and diprotons. 



in the reactions that lead to ^He. Hence essentially no 

dineutrons arc produced. The amount of additional "^He 
produced depends on the temperature at which dipro- 
ton production occurs. If the diproton is lightly bound 
the temperature will be too low for further nuclear pro- 
cessing except for the decay to ^H. If the diproton is 



tightly bound the temperature can be high enough for 
further nuclear processing to *He. In either case diproton 
production docs further reduce the hydrogen abundance. 
The second phase of ^He production occurs only if 



fu 



> 



10^ 



(G - 1.113) 



13/3 



(16) 



A typical situation is shown in figure 9. Here G = 1.3 
and fw-e = 10^^. Initially weak interactions convert 
neutrons to protons. When the Universe is 2 s old the 
temperature is 8 10^ K, which for Qd = 16 MeV is low 
enough for ^H production to occur. The ^H is quickly 
converted to *Hc, so that by age 10 s, this initial phase 
of nucleosynthesis has finished. A second phase of nu- 
cleosynthesis occurs at age 500 s, when the temperature, 
T = 6 10^ K, is low enough for production of diprotons, 
which have a binding energy of 1.8 MeV. The beta decay 
life time of the diproton is 100 s, and hence diprotons de- 
cay to ^H before significant cooling by expansion occurs. 
The temperature is sufficiently high that the ^H is con- 
verted to ^He. If the diproton binding energy was lower, 
then diproton production would occur at lower temper- 
ature and only ^H would be made in the second phase 
of nucleosynthesis. Note the small amount of neutrons 
released in the second phase of nucleosynthesis. These 
are produced by the reaction sequence d + d ^ t + p, 
followed hy t + t -^'^He + n + n. 

The dependence of the final value of Xh on G is shown 
in Figure 10 for different values of fw-e^ ranging from 
10^^ to 10^^. The value of f^^e at which Xh is signif- 
icantly reduced decreases with increasing G, due to the 
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FIG. 10: Dependence of the final hydrogen abundance on G 
for different values of fm-e ranging from 10^" at top to 10^^ 
at bottom. 



tighter binding of diproton reducing its rate of photode- 
struction. Figure 11 shows how the final value of the 
abundance depends on G for the same range of fw-e- 
Clearly a major difference from standard BBN is in the 
amount of that can be produced when the diproton 
is bound. To understand why consider the specific case 
G = 1.2. The binding energy of is then about 10 
MeV which means that it can be produced very early 
on in the big bang. Subsequent reactions reduce the 
abundance by making "^He and ''He. Most of the ''He 
is produced very quickly (90% is produced by t = 8 s). 
On the other hand, tlic binding energy of the diproton is 
relatively small, 0.7 MeV. Hence the temperature must 
drop to about 1.4 10® K before significant production 
can begin. This occurs at t = 100 s. The diproton beta- 
decay life time is about 10'^ s. Hence the diproton abun- 
dance increases during the first few thousand seconds, 
and then it decays to ^H. The temperature (< 3 10^ K) 
is now too low for further reactions involving destruction 
of ^H. Hence, depending on the value of fw-e, significant 
amounts of ^H can be produced. 

Figure 12 shows the dependence of the final ''He abun- 
dance on G for a range of fw-e values. It can be seen that 
there arc many combinations of G and fw-e for which 
complete conversion to ''He does not occur. 

For fw-e = 10^* the final H abundance is greater than 
10% of the standard value for G < 1.5. Hence signif- 
icant amounts of H remain even when the strong force 
coupling constant is 50% greater than the current value. 
In general, the final H abundance is greater than 0.075 
provided fw-e < 6 10^y{G - 1.065)^. 
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FIG. 11: Dependence of the final abundance on G for 

different values of fw-e- The thin line is for fw-c ~ 0. The 
thick lines are for fw-e ranging from 10^^ (bottom) to 10^^ 
(top). 
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FIG. 12: Dependence of the final *He abundance on G for 
different values of fw-e ranging from 10^® (bottom) to 10^^ 
(top). 



8 




1.0 1.1 1.2 1.3 1.4 1.5 



G 

FIG. 13: The solid lines are contours on which the final H 
mass fraction is 0.075 (thick) and 0.375 (thin). The thin and 
thick broken lines are the contours on which the final mass 
fraction is 0.001 with and without the pp + n and nn + p 
reactions, respectively. The dotted line is the second '*He 
production phase boundary. Below and to the left of this line 
the second phase does not occur. 

We now consider inclusion of pp + n — > d + p and 
nn + p ^ d + n as instantaneous reactions. In gen- 
eral these reactions have small effects on the final abun- 
dances, primarily because most of the neutrons have been 
depleted by the p + n reaction before the temperature 
has dropped sufficiently for diproton and dineutron pro- 
duction to occur. A small amount of neutrons are pro- 
diiccd during the second phase of nucleosynthesis by the 
p + p pp + 'J, pp d + + v, d + d t + p, 
t + t -^'^He + n + n sequence of reactions. These neu- 
trons can then react hy n+pp ^ d + p. The net result is 
a small increase in the final H abundance. Also since the 
n + pp ^ d + p reaction is assumed instantaneous, is 
produced earlier than by diproton decay alone. Provided 
the temperature is high enough, this leads to a decrease 
in the final abundance. 

Figure 13 summarizes the results of this section. The 
thicker of the solid lines is the contour on which the final 
H mass fraction is 0.075. The thinner solid line is the 
contour for final H mass fraction equal to 0.375, which is 
approximately half the standard BBN value. The broken 
lines are contours on which the final mass fraction is 
0.001. The thin solid line is the second *He production 



phase boundary. Below and to the left of this line the 
second phase does not occur. 

VI. CONCLUSIONS 

We have addressed some aspects of the effects of larger 
than standard values for the strong force coupling con- 
stant on nucleosynthesis during the hot big bang. For 
relative strong charge G > 1.065, both the diproton 
and dineutron are bound. We have estimated the beta- 
decay time scales from the ft factors for the analog nu- 
clei ^^O and ^^C. Assuming that the rate of the reac- 
tion p + p pp + J can be parameterized by multi- 
plying the rate of the reaction p + p d + + v 
by a factor fw-e-, and that the rate of the reaction 
n + n ^ nn + 7 is then related to that for the reac- 
tion pp + ^hy neglecting the Coulomb repulsion, 
we find that significant amounts of H remain provided 
fw-e < 6 10""^^/ (G — 1.065)^. By comparing similar re- 
actions, we estimate that fw-e ~ 10^^, which gives a cor- 
responding limit of G < 1.5. The primary reason for the 
survival of hydrogen is that the diproton and dineutron 
are always less tightly bound than the deuteron, which 
is a consequence of the spin-dependent part of the nu- 
clear force. Photodestruction reactions prevent buildup 
of diprotons and dineutrons before the neutrons are de- 
pleted by deuteron formation. Diprotons can be formed 
once the temperature has dropped sufficiently. These 
diprotons are converted to deuterons mainly by beta de- 
cay with possibly a contribution from the pp + n ^ d+p 
reaction. This can lead to much a larger abundance 
than in the standard BBN. 

Our main result is that the existence of bound diproton 
and dineutron nuclei does not necessarily lead to com- 
plete conversion of hydrogen to helium in the big bang. 
Instead there are parameter ranges for which significant 
amounts of hydrogen remain. We estimate for reason- 
able values of the factor by which the p -\- p ^ pp -\- ^ 
rate is enhanced relative to the p + p-^d + e^-\-v rate, 
the final hydrogen abundance is greater than 50% of the 
standard BBN value for increases in the strong force cou- 
pling constant less than about 50%. Anthropic limits on 
the strong force strength from BBN are indeed weak. 
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